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a\ 
in 

We study the flux-flow Hall effect and thermomagnetic transport near the upper 

J> ■ critical field Hc2 in extreme type-II superconductors starting from a suitable general- 
ly , 

■ ization of the time dependent Ginzburg-Landau equations. We explicitly incorporate 

■ the effects of backflow into the calculations of the local electric field and current, lead- 

<N : 

■ ing to a current which is properly divergenceless. The Hall conductivity calculated 



from this current agrees with other mean-field calculations which assume a uniform 



• applied electric field (the Schmid-Caroli-Maki solution), thereby vindicating these 

O ■ simplified treatments. We then use these results to calculate the transverse ther- 

momagnetic effects (the Ettingshausen and Nernst effects). The effects of thermal 
■ fluctuations and nonlocal elasticity of the flux lattice are incorporated using a method 

recently developed by Vecris and Pelcovits [G. Vecris and R. A. Pelcovits, Phys. Rev. 
B 44, 2767 (1991)]. Our results, taken together with those of Vecris and Pelcovits, 
provide a rather complete description of the transport properties of the flux lattice 
state near Hc2-, at least within the framework of time dependent Ginzburg-Landau 
theory. 
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I. INTRODUCTION 



The thermodynamic and transport properties of the mixed state of type-II supercon- 
ductors continue to attract the interest of theorists and experimentahsts ahke, due in large 
measure to the unusual transport properties of the high temperature superconductors. High 
transition temperatures, short coherence lengths, and large anisotropies conspire to produce 
enhanced thermal fluctuations in these materials, which can significantly modify the mean- 
field phase diagram; we refer the reader to Ref. for a detailed discussion of these effects. 
These fluctuations are also apparent in the transport properties, as they lead to a broadened 
resistive transition in the flux flow regime near Hc2 (when pinning is unimportant), and to 
thermally assisted flux flow at lower temperatures (but away from the putative vortex-glass 
transition Indeed, if we had a detailed theory of the transport properties in the presence 
of fluctuations we could in principle use this to infer properties of the equilibrium phases. 
So far, most of the theoretical work on transport properties has focused on understanding 
the behavior of the longitudinal conductivity of the flux lattice. However, it is really the 
Hall effect which represents the greatest challenge to our understanding of the dynamics of 
the vortex lattice in superconductors, as evidenced by the experimental observation that the 
Hall conductivity changes sign upon entering the mixed state in the high T^. superconductors 
a feature which is at odds with the classic theories of vortex motion in superconductors 
fSHlOl . Motivated by these observations, in this paper we re-examine the theory of the Hall 
effect in the mixed state near Hc2 using a variant of the standard time dependent Ginzburg- 
Landau (TDGL) theory. By incorporating the effects of thermal fluctuations, the nonlocal 
elasticity of the flux lattice, and backflow, we have consolidated the results of several pre- 
vious authors into a rather complete theory of the Hall effect near Hc2 (at least within the 
TDGL framework). As a byproduct, we also study transverse thermomagnetic effects such 
as the Ettingshausen effect and the Nernst effect. 

As this paper is in a sense a consolidation of the results of several different authors, it 



is appropriate to first briefly review the history of the subject. Schmid |11| first derived a 
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set of TDGL equations from the microscopic Gorkov equations. From these equations he 
was able to calculate the flux flow conductivity both near Hc2 (by solving the linearized 
equations) and near Hd (for a single vortex). The behavior near Hc2 was obtained by 
assuming that the applied electric field E was constant in space; the flux lattice is then 
effectively "boosted" by a velocity v = E x B/S^, with B the induction field. Similar 
methods were also used by Caroli and Maki |]12| to study both the dirty and clean limits. 



We will henceforth refer to this solution as the Schmid-Caroli-Maki solution. Unfortunately, 
the local current (which includes the normal current plus the supercurrent) obtained using 



this method is not divergence free, as was first pointed out by Thompson and Hu []13 
To obtain a current with zero divergence, it is necessary to incorporate backfiow currents; 
however, the backfiow has zero spatial average, so that the spatially averaged conductivity 
calculated using this method agrees with the Schmid-Caroli-Maki result. These calculations 



have recently been taken one step further by Vecris and Pelcovits [Q, who studied the 
effects of the elastic fluctuations of the flux lattice on the conductivity. Starting from the 
TDGL equations, these authors calculated the local current with backfiow, and incorporated 
the elastic fluctuations by using a dynamic generalization of of the formalism developed by 
Brandt ITHl for the static flux lattice. 



In all of these cases the TDGL equations which were employed had a purely real re- 
laxation time, and therefore exhibited a type of "particle-hole" symmetry which leads to a 
Hall conductivity which is identically zero [T^^T8|. To obtain a nonzero Hall conductivity 
one needs to generalize the TDGL equations by allowing the relaxation time to be com- 
plex. The imaginary part of the relaxation time might result from either considerations of 
Galilean invariance [l^Jl^ , or from microscopic considerations, such as fermi surface curva- 



ture Maki [|I^ and Ebisawa ||21| have used TDGL equations with a complex relaxation 
time to calculate the Hall conductivity using the Schmid-Caroli-Maki method (i.e., with- 
out backfiow). These equations have also been used to study the fluctuation Hall effect for 
temperatures T > Tc2 [ P^ , PU[ . More recently, one of us (A.T.D.) has used the generalized 



TDGL equations to study the dynamics of a single vortex (i.e., for fields close to Hd) 118 



In this paper we calculate the transport properties for the mixed state of type-II super- 
conductors starting from a set of TDGL equations which have a complex relaxation time. 
The results of this paper therefore complement the single vortex results obtained in Ref. 
||18|| . The paper is organized as follows. In Sec. II we calculate the longitudinal and Hall 
conductivities in mean-field theory by explicitly including backflow, thereby extending the 
work of Thompson and Hu and Vecris and Pelcovits. One of the important results of this 
section is that the backflow current, while important in ensuring that the total current has 
zero divergence, does not contribute to the spatially averaged Hall conductivity (which is 
of experimental relevance). Hence, the Hall conductivity which we obtain agrees with the 
result which would be obtained using the Schmid-Caroli-Maki method. We also briefly dis- 
cuss the relevance of our results to the issue of the sign change of the Hall conductivity in 
the mixed state of the high superconductors. In Sec. Ill we calculate the Ettingshausen 
and Nernst effects in the presence of backflow. Our derivation, which utilizes a recently 
discovered "virial theorem" for the equilibrium Ginzburg-Landau equations, is exact within 
mean-field theory. The effects of elastic fluctuations of the flux lattice are considered in Sec. 
IV, which follows the work of Vecris and Pelcovits. We find that the amplitude fluctuations 
in the flux lattice phase suppress the flux flow conductivities below their mean-field values. 
Nonlocal effects are extremely important in setting the scale for these fluctuations. 

II. MEAN FIELD THEORY 

The TDGL equations consist of an equation of motion for the order parameter ip, 

ill + tl2){dt + t^)^ = {-- iA.f^lj + ^ - I^P^, (2.1) 

K 

along with Ampere's Law, 

Vxh = J, (2.2) 

where h = V x A is the local magnetic induction field. For the current we adopt a two fluid 
model, so that J = J„ + J^. The normal current is 
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■ E, (2.3) 



where the electric field is expressed in terms of the potentials as 

E = --V$-9tA. (2.4) 

In the normal current we include both the longitudinal and the transverse response of the 
normal carriers; the "■" which appears in Eq. ( p.3|) indicates a tensor product, with a^'^^ the 
normal state conductivity tensor, 

, 7,1. (2.5) 

\ xy ^ XX . 



The signs used in cr^"-' are appropriate for positive carriers. The supercurrent is given by 

-ipVifj*) -l^l'^A. (2.6) 



xy 

1 



2Ki 

These equations are written in dimensionless variables such that lengths are scaled by the 
magnetic penetration depth A, time is scaled by h/2mC,'^ with C, the coherence length, mag- 
netic fields are scaled by V2Hc with He the thermodynamic critical field; k = X/^ is the 
usual Ginzburg-Landau parameter. As an aid to the reader, important results will be ex- 
plicitly expressed in both dimensionless and conventional units. The quantities 71 and 72 
are the the real and imaginary parts of the dimensionless order parameter relaxation time. 
The scalar potential is denoted by $; the difference between the scalar potential and the 



electrochemical potential will be ignored here (see Refs. and \T^] for a more extended 
discussion). Since in equilibrium we will assume that we have local charge neutrality, out 
of equilibrium any charge density must be 0{v), with v the velocity of the vortex lattice; 
therefore the time variation of the charge density is O(f^), and will be neglected in the 
spirit of the linear response calculation of this paper. As a result, the total current must be 
divergenceless; i.e., V ■ (J„ + Js) = 0. 

Before attempting to solve the TDGL equations, it is useful to first simplify them some- 
what. To do this, we write the order parameter in terms of an amplitude and a phase. 
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ipir, t) = /(r, t) exp[i(p{r, t)]. In terms of the gauge invariant quantities Q = A — Vip/ k and 
P = $ + dt(p, the magnetic and electric fields are 



h = VxQ, (2.7) 



E = --VP-dtQ, (2.8) 



and the supercurrent is 



-fQ. (2.9) 



The real part of Eq. ( |2.1| ) is 



lidj - l2Pf = ^VV -Q'f + f- f, (2.10) 



while the imaginary part is 



l2dtf + iiPf + -/V ■ Q + -Q ■ V/ = 0, (2.11) 

K K 



and Eqs. (Pg)-(P;B) become 



V X V X Q = ■ (--VP - dtQ) - fQ. (2.12) 

An explicit equation for P may be obtained as follows. First, multiply Eq. (|2.11|) by /; the 
gradient terms can be combined as V ■ J^; then use the fact that V ■ = —V • Jn- We 
finally obtain 

-V ■ ■ (--VP - dtQ)] + 7if P + l2fdtf = 0. (2.13) 

We begin by calculating the local electric field for the moving fiux lattice. First, we 
assume that the lattice translates uniformly, so that /, Q, and P are only functions of 
r — vt. Therefore we replace all time derivatives in Eqs. ( |2.1CI| ), (|2.12| ), and ( |2.13| ) by — v ■ V. 
Second, as we are concerned with linear response in this paper, we keep only terms of order 
the fiux lattice velocity v. In this spirit, we expand all quantities in powers of the velocity. 



with the order of expansion denoted by a superscript: / = /^°-' + /*-^\ Q = Q^^-' + Q^^-', where 
/(^) and Q(i) are 0{v). Note that P is 0{v), since the electric field vanishes in equilibrium. 
The 0(1) equations are simply the equilibrium Ginzburg-Landau equations. The electric 
field can therefore be written as 



E 



-V X 



1 



h(o) _ V -P - V ■ Q 



(0) 



K 



(2.14) 



Upon averaging over the volume V of the sample, we find for the spatially averaged electric 
field 

1 



(E) 



V 



d'^r E(r) = -v x B, 



(2.15) 



since the average of the gradient term in Eq. ( p.l4|) can be converted to a surface term 
which vanishes at the boundaries; B = (h^^-') is the (equilibrium) macroscopic induction 
field. Although the gradient term in Eq. (|2.14| ) does not contribute to the spatially averaged 
electric field, it does contribute to the local electric field. We therefore need to calculate 
P/ K — \ ■ Q'-^-'; an equation for this quantity follows from Eq. ( |2.13| ): 



V 



^(") . V ( -P - V ■ Q(°) 
= -V 



(0) (ip_v.Q(o) 

K 



(")-(vxhW) + ^a;Wv 



+ kVw^°^v ■ Q 



)(0) 



(2.16) 



where for simplicity we have introduced tu^^-' = {f^'^^Y- The last term on the right hand 
side of Eq. ( |2.1(j| ) can be further simplified by noting that from the equilibrium equations 
we have v ■ (tu^'^^Q,^^^) = V ■ (v x h^^-'). Simplifying the derivatives on the left hand side, we 
finally arrive at 



ai^) (ip - V . Q(°)) - . 



2^^^(0) ( ip _ V . Q 



(0) 



Vj, 



(2.17) 



where we have defined 



j(r) = (t(") ■ [v X 



h(o)(r)] - «:2^iv X h(0)(r) + '^J°){r)^r. 



(2.18) 



Next, define the local deviations from the average equilibrium values of the magnetic induc- 
tion field and the square of the order parameter as 



5h{r) = h(°)(r) -B, 



(2.19) 



6uj = uj^''^-{J'^), 



(2.20) 



so that ((5h(r)) = {6uj) = 0. When these expressions are substituted into Eq. (|2.18| ), there 
will be a constant piece which can be discarded as it will not contribute to Eq. ( |2.17| ). Noting 
that for the equilibrium state 6h = —{6uj/2k)z, which is correct to 0{6uj'^) we then see 
that it is possible to write Eq. (p.l8|) in the following form: 



j(r) = a^") ■ [v X 6h{r)] - —6u{r){z x v) + —6u{r)v 



7iK 



72 



(2.21) 



Eqs. ( p.l4|) , (|2.17| ) and ( p.21|) will together determine the local electric field, and therefore 
the local normal current. 

The solution to Eq. (|2T7D is 



-P(r) - vQ(°)(r)= - / rfVG(r,r')V'-j(r') 
= J dVj(r') ■ V'G'(r, r'), 



(2.22) 



where G{r, r') is the Green's function which satisfies 



ai^V - Ai^^°Hr)l G(r, r') = ^(^^(r - r'), 



(2.23) 



and where in the second line of Eq. ( 2.22 ) we have integrated by parts and neglected a 
surface contribution. With this solution it is possible to calculate the normal current. First, 
take a gradient of Eq. ( p.23|) and multiply by the normal state conductivity tensor; after 
using several vector identities, we find 

a(-) . V - V • Qo) = I dV j(r')a(^^V • V'G'(r, r') 

-V X y rfV'j(r') X ■ V'G(r,r')] . (2.24) 

The flux lattice state is not translationally invariant, so that G{r, r') is not a function of the 
coordinate difference alone. However, sufficiently close to ifc2 the order parameter amplitude 
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is small, and we can replace uj^^\r) in Eq. (|2.23|) by its spatial average {ou^^^), which is 
correct to 0{6uj). Within this approximation, the Green's function ^(r, r') = G{r — r'); 
then V • V'G'(r, r') = —V^Glr — r'). Combining this result with Eq. (|2.23|) for the Green's 
function, we find that Eq. (|2.24|) becomes 



. V - V ■ Q(°)) = -j(r) - / rfVG(r - r')j(r' 



V X / d\'i{r') X ■ V'G(r - r' 



(2.25) 



The second term on the right hand side of Eq. ( p.25| ) will generally be quite small near the 
transition, as j itself is 0{6h); this is then multiplied by (ujq), rendering the second term 
doubly small near the transition. We will therefore drop this term in what follows. The 
normal current is obtained by combining Eq. ( p.25|) with the definition of the normal current, 
Eq. ( |2.3|) , along with the expression for the electric field, Eq. ( p.l4|) ; the final result is 



Jn(r) 



-a(") ■ [v X B] - ^6uj{r){z x v) + '^6uj{r)v 
+V X / c/V j(r') X ■ V'G(r - r' 



(2.26) 



As a check on our result, we note that the spatial average of the last three terms on the right 
hand side of Eq. ( |2.26| ) is zero, so that (J„) = — cr^"^ ■ (v x B) = a^"'^ ■ (E), as required. It is 
also straightforward to show that when a^^^ =72 = 0, Eq. ( p.26| ) reduces to the analogous 



results derived by Thompson and Hu |]T3[ and Vecris and Pelcovits |Tj] in two and three 
dimensions, respectively. 

We next calculate the linearized supercurrent for the moving flux lattice. This calculation 
is most conveniently carried out in the symmetric gauge, rather than the Landau gauge used 
by Thompson and Hu [|l^]. Following Vecris and Pelcovits |T^, we start by using a postulated 
solution for the order parameter of a uniformly translating flux lattice, 

Bk 



^piir{t)) = exp 



4 



g{x{t) + iy{t)), 



(2.27) 



where the subscript / indicates a linearized solution, r(t) = r — vt is the coordinate in the 
moving frame, and g{x{t) + iy{t)) is an analytic function (to be specified later). In the 



presence of an electric field there are corrections to k of ©(t;^), which we drop in the spirit of 
our linear response calculation |[T^. Substituting Eq. ( |2.27] ) into the first TDGL equation, 
Eq. and dropping terms proportional to 'Wg/ g , we find the required symmetric gauge 

potential near Hc2, 



A = —z X r 
2 



z X V H V. 

2 2 



(2.28) 



We use an analytic function, g{x{t) + iy{t)) |]T1|,|T3|, appropriate for a translating flux-line 
lattice with flux lines located at and parallel to the z-axis. 



N 



g{x{t) + iy{t)) = \[{x{t)-x, + i[y{t)-y,]}. 



(2.29) 



This form of the order parameter has zeros at the instantaneous vortex positions r,^. We 
therefore have for the square of the amplitude of the order parameter. 



2 ^ ' 



uJi{r{t)) = exp 
and for the phase of the order parameter 

N 

(^(r(t)) = ^tan-i 



N 



n |rW 



(2.30) 



y{t) - yu 

X{t) - Xy 



(2.31) 



For the gauge invariant vector potential we then have 



Q= A Vy? 



B ^ 7iK „ 72K 1 Jl, i X [r(t) 
-2; X r —z X V H V 2^ 



2 2 
z X VcJi 



2KCi;/ 

The linearized supercurrent 

J. 



2; X V H V. 

2 2 



(2.32) 



-cuQ is therefore 



+ ^cK^ X V) 



72 



-CL^iV. 



(2.33) 



where J^^-* = i x Vuji/2k is the uniformly translating equilibrium supercurrent. 

The total current J = J„ + Js is obtained by adding our expression for the normal 
current, Eq. ( |2.26| ), to our expression for the supercurrent, Eq. ( ^.33| ), which we obtained 
above. Using v = (E) x B/_B^, we find 
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J(r) = a • (E) + J(°) + V X J ciV j(r') x [a^") ■ VG{r - r')] , (2.34) 
where the conductivity tensor is 

\ '^xy ~r 2_B ^ xx ~ 2B ' 

The spatial averages of the last two terms on the right hand side of Eq. ( |2.34] ) are zero, so 
that (J) = o ■ (E); we also have V ■ J = 0, since V ■ J^"-* = 0. We have therefore found a 
current which is properly divergenceless, as required. The various terms on the right hand 
side of Eq. (|2.34|) also have simple interpretations — the first term is a uniform transport 
current, the second is the uniformly translating equilibrium supercurrent, and the last term 
is the backflow current. This form of the local current was first obtained by Thompson and 
Hu ||13| for the case cr^!^^ = 72 = 0; our result is a generalization to the situation in which 
there is particle-hole asymmetry. 

In mean-field theory the Abrikosov value for (cu^'^)) is (see Ref. [Q, for instance) 



m Eri - B 



where Pa = ii^^^^)'^) which is 1.16 in mean-field theory for a triangular flux lattice 
This leads to the following expressions for the conductivities in mean-field theory, in 



conventional units: 



in) , 71"^ 1 Hc2-B 



_ In) I /I"" ^ -"-'c^ fry oyN 



and 



Notice that the conductivities have contributions from both the normal carriers and from the 
vortex motion. The real part of the order parameter relaxation time, 71, is always positive, 
so that this contribution is additive for the longitudinal conductivity. However, the sign of 
72 is most likely determined by microscopic considerations if 72 > 0, then it is possible 
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for the Hall conductivity to change sign in the mixed state. Further microscopic calculations 
are needed to determined if this is the source of the sign change which has been observed in 
the high Tc superconductors. 

It is also possible to calculate the corrections to the local magnetic field for a moving 
flux lattice. We can do this by expressing the local current as a curl: 



V X h(r)= J(r) 

= -V X [z (Jt X z) ■ r] - V X 



zJ'^\r)/2K 



+V X y dV j(r') X [(t(") • VG{r - r')] , (2.39) 

where Jj = cr ■ (E) is the uniform transport current. Integrating, we obtain [h(r) = h{r)z]: 

h{r) = /i(°)(r) - (Jt X z) ■ r + J rfV j(r') x [ct^") • VG(t - r')] , (2.40) 

where h^^^ = B — 5iu/2K is the equilbrium local magnetic field. This is a generalization 
to the particle-hole asymmetric case of the result of Vecris and Pelcovits (see Eq. (2.23) of 
Ref. [n|). As noted by these authors, the second term in Eq. ( p^.40| ), which grows linearly 
with distance within the sample, is typical of magnetostatics problems in the presence of a 
uniform current density. 

To summarize our results so far, we have explicitly calculated the total current for a 
moving flux lattice starting from the generalized TDGL equations. This current has zero 
divergence; however, the conductivities calculated from this current are identical to those 
which would be obtained by using the Schmid-Caroli-Maki solution. 



III. THERMAL TRANSPORT 



The moving flux lattice not only produces dissipation but also transports energy, in a 
direction parallel to its velocity. In order to calculate the transported energy in the mixed 



state, we start from the expression for the energy current which is due to Schmid [pJ|] : 

'V 



2Exh-2ExB + 2 



- iA 4j(dt - i^W + c.c. 
n J 



(3.1) 
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where J'* is in units of (-ff^/47r)(^/2m)(/t^/A) (i.e., units of energy per unit volume times 
velocity). The second term, which was not considered by Schmid, is necessary in order 



to subtract out the contribution from the uniform background field B = (h) ||23|. This 
expression may be greatly simplified by using a sequence of transformations which were used 
in Ref. |T^; these are reproduced here for completeness. First, in terms of the potentials P 
and Q, along with the order parameter amplitude /, to 0{v) Eq. ( p. I]) becomes 



1 




— VP + V 






2 




+- 


[-i(v 


K 





^(0)^ 



V V/W)(V/ 



q(o)) _ 2E X B 

pQ(0)(j(0))2l 



(3.2) 



Using V X V X Q(°) + (/(°))^Q(°) = 0, the first and last terms in Eq. (|3^ ) may be combined: 
(VP) X (V X Q(°)) + PV X V X Q(°) = V X (PV x Q(°)), (3.3) 



where a vector identity has been used. The second term on the left hand side of Eq. ( p.2|) 
may be written as 

(v ■ VQ(°)) X (V X Q(o)) = V X [(v ■ Q(°))V X Q^^)] + v(V x Q^o))^ - (v ■ Q(°))V x V x Q(°), 

(3.4) 

where we have again used several vector identities. Combining Eqs. ( |3.2| )-(p^, we have 



r= 2v X 

+2 



1 



2E X B 



(v ■ V/(°))(V/(°)) + (/(°))'(v ■ QW)QW + v(/iW)^ 



(3.5) 



Therefore the backflow terms only appear in the first two terms on the right hand side of 
Eq. ( |3.5| ). However, when we calculate the spatially averaged heat current, the first term on 
the right hand side of Eq. ( p. 51 ) can be converted into a surface term, which vanishes; the 
second term yields —2(E) x B = — 2P^v; therefore the backflow corrections do not enter 
into the calculation of the spatially averaged energy current. After performing the spatial 
average, we find that (J^) = nU^w, where n = k,B/2tt is the vortex density (n = P/0o iii 
conventional units), and where is the transport energy per vortex; we have 
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1 (V/(o))2 + (/(o))2(g(o))2^2(/i(°))2 



2B\ (3.6) 



(a factor of 1/2 arises from an angular average). The first two terms are half of the kinetic 
energy of the superfluid, while the third term is the magnetic field energy. Recently, Doria 
et al. have proved a "virial theorem" for the equilibrium Ginzburg-Landau equations ||24|| , 
which shows that the integral which appears in Eq. (|3.6| ) is precisely equal to 2H ■ B. 
Therefore we find that 

nU^ = -2(B - H) ■ B = -SttM ■ B, (3.7) 

with M = 47r(B — H) the spatially averaged equilibrium magnetization of the sample. 
Therefore, the transport energy in dimensionless units is 

Us= A-kM 

K 

An K- H 

while in conventional units we have 



(3.9) 



AtiI3a(2k'^ - 1) 



where we have used the mean-field expressions for the magnetization [^. Using the lin- 



earized microscopic theory near Hc2-, Maki obtained the result = —(f)oMLD{t), where 



t is the reduced temperature; L£){t) ^ 1 in the dirty limit near Hc2, so that our results agree 
in this limit. Note that our result is much more general, as the derivation did not invoke the 
assumption of linearity of the order parameter, but only the assumption of linear response 
in the flux lattice velocity. Therefore our result holds for the entire mixed state, and not 
just near Hc2- 

The thermomagnetic transport coefficients for a superconductor in the mixed state are 
discussed in Ref. [|TB|, to which we refer the reader for details. The Nernst coefficient is 
defined as u = Ey/H{dT/dx), under the conditions of Jx = Jy = dT/dy = 0. Introducing 
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the transport coefficient axy through (J'') = axy(E) x z, then it is possible to write the 
Nernst coefficient a.s u ^ {l/TH){axy/ (Jxx)-, where a^x is the full conductivity (including 
both the normal state and flux flow contributions); but from the above discussion we see 
that Uxy = U^/(f)Q. We therefore flnd for the Nernst coefficient, 

(poTH axx 

The Ettingshausen coefficient is deflned as £^ = {dT/dy)/HJx under the conditions Jy = 
Jy = dT/dx = 0. Using the Onsager relations, it is possible to show that S = Tu/kxx, 
where the thermal conductivity. 

IV. FLUCTUATIONS 

Having calculated the transport properties in the mean-fleld regime, we now turn to the 
effect of thermal fluctuations of the flux lattice on the transport properties. To do this, 
one flrst assumes that the flux lines are located at rj,(z) = R^, + s^{z), where {Ru} are the 
positions of the flux lines in mean-fleld theory (which form a triangular lattice), and {s^{z)} 
are the deviations from the mean-fleld positions. Expanding about the mean-fleld solution, 
and then taking the continuum limit by replacing Si^(z) by s(r), the free energy becomes 
F = Fq + Ff^i, where Fq is the free energy of the mean-fleld Abrikosov state and F^i is the 
elastic free energy given by 

^ei = lJ s,(-k){[cn(k) - cM]hkj + 5ij[cMkl + C44(k)A;2]}s,(k), (4.1) 

where cn, C44, and cee are the uniaxial compression modulus, tilt modulus, and shear modu- 
lus, respectively. The derivation of the nonlocal elastic moduli from Ginzburg-Landau theory 
was flrst carried out by Brandt |T5[ ; his results have recently been generalized to the case of 



anisotropic superconductors by Houghton, Pelcovits, and Sudb0 . The current must now 
be averaged with respect to an ensemble specifled by the elastic free energy; as shown by 



Vecris and Pelcovits WM, this is equivalent to replacing the spatial average of the square of 
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the mean-field order parameter, {uj{r)), which appears in the expression for the mean-field 
conductivities, Eq. ( |2.35| ), by the spatial and ensemble average of the square of the order 
parameter, ((a;(r))th), where (■ ■ ■)th is the ensemble average (we will drop the superscript 
on uj for simplicity). Therefore, in the presence of thermal fluctuations the longitudinal 
conductivity (first obtained by Vecris and Pelcovits) becomes, in conventional units, 

^.x = aW+7ie*((^(r))th)/5, (4.2) 

and the Hall conductivity becomes 

^.. = <^-72e*((o.(r))th)/i?. (4.3) 

A generalized form for the square of the order parameter in the mixed state was first 



suggested by Brandt [T^ for isotropic superconductors. Here we will use an anisotropic 
generalization, appropriate for superconductors with an anisotropic effective mass tensor; 
this generalization is straightforward and is implicit in the work of Houghton, Pelcovits, and 



Sudb0 The order parameter is 



f r r dk r d'^k, pjk [r-r^(z')] ~l 

^ iv,B) exp {-4. E / / 1 L wr^iiTiwn^} • (^-^i 

where N{B) is a magnetic field dependent normalization constant, 7^ = m/mz is the mass 
anisotropy, with m the effective mass in the plane and the effective mass along the z-axis, 
k_L = {kx, ky), BZ denotes an integration over the first Brillouin zone, k"^ = 2(1 — &)/^af) 
conventional units, with the in-plane coherence length, and b = B/Hc2 is a dimensionless 
magnetic field. As argued by Brandt, this form of the square of the order parameter has the 
proper second-order zeros at the vortex positions {r,,}, and reduces to the correct forms for 
both large and small inductions. Expanding the exponent to first order in 8,^, and taking 
the continuum limit, we obtain 

dk, f d^k^ k ■ s(k) e-^^"- 



iu[r) = uJo[r) exp 



271 Jbz (27r)2 kl + -i^kl + kl 



(4.5) 



where ujq{y) is the square of the mean-field order parameter with zeros at {Rj/}, and n 
is the vortex density; the fc^ integration is now over a circular Brillouin zone of radius 
kBz = (47m) = (26)^/^/^ab. Performing the thermal average, we obtain 
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(u;(r))th = uJo{r)e , 



(4.6) 



where W is a suppression factor, given by 



/bz (27r)2 (kl + -i^kl + klf 
The fluctuation propagator is given by 



fcifcj(si(k)sj(-k))th. 



(4.7) 



(5,(k)s,(-k))th = kBT[PmT(y) + P'^Gi(k)] 



(4.8) 



where P^^ = Sij — kikj/kj_ and Plj = kikj/kj_ are the transverse and longitudinal projection 
operators, respectively, and where the transverse and longitudinal propagators are given by 

1 „ . 1 



C66(k)A;]_ + C44(k)A;2' Cii{k)kl + Cu(k)kl' 

Carrying out the implicit summation in Eq. ([4.7|) , we obtain 

d'^k , A;2 1 



(4.9) 



Z J Zn Jb: 



/bz (27r)2 {kl + ^^kl + klY Cn(k)A;i + Ci^{k)kl 
The spatial average is trivial; our final result is 

((u;(r)),h) = (u;o(r))e-^. 



(4.10) 



(4.11) 



with {ujq{y)) the spatial average of the square of the mean-field order parameter given in 
Eq. (^. 

In order to calculate the supression factor W we first rescale the momenta by the Brillouin 
zone radius by introducing a dimensionless variable q = k/k-Qz- Then we have 
1 1 ^2 ^oo .1 I 



W 



/ dq^ / dqj_ 
Jo Jo 



9 ' 



(4.12) 



TT At/cbz 47r Jo """^ Jo '^^ [qi + 72^2 + m2]2 cii(q)g5. + C44(q)g2 
where = ^Q/lOvr^fc^T is a thermal length |jl|, and where m| = {k^/k-QzY = (1 ~ ^)/^- 



The nonlocal, anisotropic elastic coefficients are given by p6|j27 

2 



ciiiqj 



/ kh 



Att \ k 



BZ, 



(g2 + m1){qj_ + 7^^^ _|_ ■y'^ml) q\ + 7^5^ + m|j 



(4.13) 
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5' / kh 



2 



1+ 



(4.14) 



with A;^ = (1 — 6)/A^, and m| = {kh/k-BzY = (1 — 6)/26k^. The integral is rather formidable, 
and we have not succeeded in evaluating it in a general form. However, progress is possible 
if we consider the limit m\ <^ 1, i.e., 1/{2h?) ^ 6/(1 — 6), which is easily satisfied for most 
of the mixed state in the high superconductors. In this limit it is possible to take to 
nonlocal limit of the elastic coefficients, i.e., set m\ = in the expressions for cn and C44, 
Eqs. (|4.13| ) and (|4.14|) . In this limit the anisotropy factor 7 only enters when multiplied by 
Qzi so it may be scaled out of the integral. The integral on can be performed analytically, 
although the result is quite complicated. The remaining integral can then be performed 
in the limits m| ^ 1 or ^ 1. We will consider these cases in turn. 

(i) m| ^ 1; i.e., 1/(2k^) ^ 6/(1 — 6) ^ 1 . In this limit, which applies to most of the 
mixed state, we obtain 



m J Ar^BZ V k^ 
v^/27rem,y/2 ^ ^5/2 



V m ; (1 -t)V2 (1 _5)3' ^^•^^) 

where we have used ^ = ^ab(0)(l— t)^^/^ with t = T/T^] e is the Ginzburg criterion parameter, 
given by 

_ 1 ( ^'Um ' _ l6n^>^\kBT,? ..... 

The Ginzburg parameter measures the relative strength of thermal fluctuations; the high 
Tc superconductors typically have (em^/m)^/^ of 0(1) or larger, an enhancement of several 
orders of magnitude over the low Tc superconductors ||l|,|26|. These fluctuations lead to 
a reduction of the conductivity and a significant rounding of the resistive transition in a 
magnetic field. 

(ii) m| <^ 1; i.e., 1/(2/?^) ^ 1 ^ 6/(1 — 6) . This is the appropriate limit for fields near 
the mean- field upper critical field ifc2- In this limit we obtain 
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\m J ArkBz \ J 
' > (4.17) 



2 V m y (1 -t)V2 (1-5)3/2 • 

Again, we note the appearance of the Ginzburg parameter e. We would also like to point 
out that W exhibits an interesting scahng behavior in this hmit; i.e., 



W{B,T) = T 



1-6) 



(4.18) 



where the scahng function J-'{x) = 0.29 x^^^^. Therefore the conductivities near Hc2 exhibit 
a scahng behavior quite simihar to the scahng behavior which is inherent in transport calcu- 
lations which use the lowest Landau-level Hartree approximation [|T6| or extensions thereof 
PE| . This scaling behavior has recently been observed in measurements on the high 
superconductors |]29[| . 



The importance of nonlocal effects in setting the scale for thermal fluctuations which 
may melt the flux lattice was first pointed out by Brandt |30|. This can also be observed in 



the amplitude fluctuations. In the isotropic limit (7 = 1) the suppression factor in the local 
limit is easily calculated; we find 

12 At (1-6)2' ("■^■'i 
which lacks the important factor of which appears in the nonlocal expression. 

There are several features of our result for the thermally averaged order parameter which 
are noteworthy. First, the fluctuations suppress the conductivity below the mean-field value; 
this should be contrasted with the behavior above the mean-field Hc2, where fluctuations 
enhance the conductivity above the normal state value. Second, notice that the integral 
in Eq. ( ^.10| ) is infrared convergent; there are no divergences in the amplitude fluctuations 



which we consider here pT|. This is in contrast to the phase fluctuations, which diverge 



with the system size [p2| , p3| , p7| . This divergence is taken as an indication of the absence of 
off-diagonal long-range order (ODLRO) in the flux-lattice state [^. The amplitude fluc- 
tuations, while suppressing the conductivity below the mean-field value, do not drive the 
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flux-flow contribution to the conductivity to zero. We therefore have an enhanced conduc- 
tivity even in the absence of ODLRO. Third, the amphtude fluctuations are longitudinal, 
unlike the phase fluctuations which are transverse. As a result, our expression for W does 
not involve the shear modulus cqq. It would then appear that the conductivities are rela- 
tively insensitive to a vortex lattice melting transition |I| J2^ , |5^ , at which the shear modulus 
would be abruptly driven to zero in crossing the liquid-solid phase boundary |3^. However, 



this observation may be signiflcantly modified once we account for vortex pinning |ll|, |36| , |37 



V. CONCLUSIONS 

To summarize, we have calculated the transport coefficients in the mixed state using a 
generalized TDGL theory. Our calculations have explicitly incorporated "backflow" effects, 
yielding a current which is properly divergenceless. However, the the results which we obtain 
are wholly equivalent to the Schmid-Caroli-Maki solution of the TDGL equations, since 
the backflow current has zero spatial average. Therefore, at least within the framework of 
TDGL theory, the backflow currents associated with vortex motion have little bearing on the 
question of the sign change of the Hall conductivity. We also calculated the thermomagnetic 
transport properties in mean-field theory, and found that under quite general circumstances 
the transport energy is proportional to the equilibrium magnetization. Finally, we find that 
elastic fluctuations of the vortex lattice tend to suppress the conductivities. Hartree-type 
approximations, which extrapolate from the high temperature phase, generally predict an 
increase in the conductivities due to fluctuations |[T6|j28|l . How to reconcile these two quite 
different approaches remains an open problem. 
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